蝎 虎 天 体 PKS 07354178 的 光 变 特性 分 析 


RE, he, ELT, FRB 
(云南 师范 大 学 物理 与 电子 信息 学 院 ， 云 南 昆明 650500) 


摘要 : 在 收集 大 量 数 据 的 基础 上 ， 用 时 间 补 偿 离 散 传 里 叶 变 换 、Jurkevich 方法 和 离散 相关 
分 析 法 分 析 了 PKS 0735+178 的 B 波 段 和 V 波 段 光 变 周 期 ,发 现 该 天 体 具 有 (4.33 士 0.41) 


年 的 光 变 周期 ， 其 中 心 黑洞 质量 的 下 限 为 0.22 x10sM 
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BL Lac 天 体 是 活动 星系 核 中 一 个 重要 的 子 类 ，BL Lac 天 体 在 观测 上 表现 出 高 光度 、 
偏振 、 快 速 光 变 以 及 非 热 辐射 等 特征 0 2，BL Lac 天 体 具有 长 周期 光 变 和 短 时 标 光 变 ， 
过 观测 和 研究 BL Lac 天 体 不 同 的 光 变 时 标 能 获得 天 体 的 中 心 黑 洞 质量 、 辐射 区 域 及 内 部 
结构 参数 等 “"。 因 此 ， 研 究 该 天 体 的 光 变 周期 非常 重要 ，。 

BL Lac 天 体 PKS 07354178 是 一 颗 红 移 为 0.424 KRZR”, 它 具 有 平 射电 谱 和 射电 爆 、 
剧烈 的 光 变 、 高 偏振 及 超 光 速 运动 等 性 质 ”"， 不 发 出 或 仅 发 出 微弱 和 间断 的 发 射线 . 众 所 
周知 ， 分 析 PKS 0735+178 光 变 周期 的 方法 很 多 ， 但 是 有 些 方 法 误差 极 大 ， 对 观测 数据 的 连 
续 性 有 极 高 的 要 求 , 并 不 适合 周期 性 研究 , 而 BL Lac 型 天 体 由 于 观测 特征 表明 无 发 射线 "”， 
要 通过 光谱 观测 获取 光谱 线 ， 并 用 其 研究 中 心 黑 洞 质量 等 内 部 结构 参数 是 不 可 能 的 ~”。 本 
文 主要 使 用 时 间 补 偿 离散 傅 里 叶 变 换 、 离 散 相 关 分 析 法 和 Jurkevich3 种 方法 对 B 波段 和 V 
波段 的 光 变 周期 进行 研究 , 并 对 比分 析 这 3 种 方法 ,其 中 时 间 补 偿 离散 傅 里 时 变换 和 离散 相 
关 分 析 法 是 初次 用 于 研究 PKS 0735+178 的 光 变 周期 ， 这 两 种 方法 对 观测 数据 的 连续 性 要 求 
低 、 结 果 准 确 ， 而 Jurkevich 方法 在 文 [9] 中 用 过 ， 并 且 它 要 求 观测 数据 的 时 间 序 列 长 、 连 
续 性 强 。 


fai m 


1 样本 和 光 变 曲线 


本 文 研究 的 数据 从 文 [11-18] 获 取 ， 收 集 从 1970 年 到 2002 年 期 间 PKS 0735+178 光学 也 
波段 和 V 波段 的 观测 数据 。 图 1 和 图 2 分 别 为 B 波段 和 Y 波段 的 光 变 曲线 。 从 光 变 曲线 可 以 
看 出 : PKS 0735+178 在 光学 波段 活动 非常 激烈 ， 在 近 30 年 的 观测 中 ，B、V 波段 的 最 大 变化 
ir 3. 5 个 星 等 ， 由 于 受 观测 条 件 的 限制 ， 光 变 曲 线 数据 不 连续 ，V 波段 将 近 9 年 没有 数据 ， 
使 周期 性 分 析 受 到 限制 。 
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图 1 PKS 0735+178 在 B 波段 的 光 变 曲线 图 2 PKS 0735+178 4E V 波段 的 光 变 曲线 
Fig. 1 The light curve of quasar PKS 0735+178 in Fig. 2 The light curve of quasar PKS 0735+178 in 
B band V band 


2 周期 分 析 


2.1 时 间 补 偿 离散 傅 里 时 变换 分 析 PKS 0735+178 的 光 变 周期 

时 间 补 偿 离散 傅 里 时 变换 方法 是 计算 光 变 周期 最 常用 的 方法 之 一 , 文 [19] 用 该 方法 分 析 PKS 
1510-089 红外 光 变 周期 , 通过 对 1 sinoat, cos œt fE Gram-Schmidt 正 交 化 ， 得 到 3 个 正 
交 向 量 ， 将 数据 投影 到 3 个 正 交 向 量 上 就 得 到 了 频谱 : 


A, =1 ; (1) 
H,-cosot | (2) 
H, =sin ot (3) 
正 交 化 后 : 
hy = aH, (4) 
h, = aH, - ayay(h,. Hi) (5) 
h, -ajH, —a,a,(h,H,)—a,a,(h,,H,) l (6) 


对 于 均匀 采样 的 数据 ， 这 个 过 程 对 应 于 用 三 维 “ 正 弦 + 常 数 ” 模 型 进行 曲线 拟 合 。 如 果 周 期 
大 于 采样 时 间 ， 且 时 间 序 列 足 够 长 ， 履 盖 所 有 的 相位 ， 则 有 


1 N 
(H,,H,) = (Hy, Hy) 854977 
^ qu 
功率 谱 强度 为 : 
UT I9) -2ajl|F(»|- c c E 


其 中 ， 下 (@ ) 为 离散 傅立叶 变换 的 功率 谱 。 在 不 均匀 采样 的 情况 下 ， 下 (@) 是 加 权 的 时 间 补偿 
离散 傅立叶 变换 : 
F(o) 2 Cf. h, +ih,)/a,v2 9) 
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在 许多 蝎 虎 天 体 的 观测 中 ， 观 测 数据 f (7, ) 的 精度 各 不 相同 。 考 虑 到 这 个 问题 ， 引 进 权重 方 


FEO; = Wy, ) 重 新 定义 内 积 ; 


(81,85) =>, ,8;(t,) g5(t;) (10) 
频率 9 处 的 强度 由 下 式 给 出 : 


I(@) 2c; +c (11) 
由 线性 回归 理论 可 知 , 0< 了 (@) < Q. JU Q - (f, f) - Yo, f) HAHEN, bit 


标准 化 因子 : 统计 量 S(@) = 1(0)/ Q. 称 这 一 量 为 谱 相关 系数 ， 对 于 所 有 频率 9 ， 


0<S(o)<1. 
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FREQUENCY (1/day) FREQUENCY (1/day) 
图 3 DCDFT 法 分 析 PKS 07354178 的 B 波段 周期 图 4 DCDFT 法 分 析 PKS 0735+178 的 Y 波段 周期 
Fig.3 DCDFT method analysis the variability Fig.4 DCDFT method analysis the variability 
period of PKS 0735+178 in B band period of PKS 0735+178 in V band 


应 用 上 述 方法 获得 PKS 0735+178 的 B 波 段 和 V 波段 周期 图 ， 由 图 3、 图 4 可 知 ，B 波段 
存在 0.62 4E. 1.01 E. 1.67 年、4.72 年 的 光 变 周期 信息 ，V 波段 存在 0. 58 年 、1. 07 年 、 
1.67 年 4.72 年 的 光 变 周期 信息 , 按 文 [19] 中 的 判 据 ，B 波段 存在 1.01 年 和 V 波段 存在 1. 07 
年 光 变 周期 信息 的 可 靠 性 最 大 ， 而 B 和 YY 波段 存在 4, 72 年 次 之 . 但 考虑 到 平均 时 标 为 


(1.11 土 0.01) 年 的 周期 很 可 能 是 由 于 地 球 绕 太 阳 所 致 , 则 光 变 周期 B 波段 为 4. 72 年 , V 波段 


Wy 4, 12 5&. 
2.2 Jurkevich 分 析 周 期 
Jurkevich 基于 对 天 文 测量 中 的 非 均 匀 测 量 问题 提出 的 一 种 统计 法 ”， 如果 观测 样本 数 


据 为 N 个 ， 久 ,为 单 次 测量 值 ， 及 为 所 有 测量 值 的 平均 值 ，V” 为 测量 数据 的 方差 ， 5? 为 测 
量 数 据 样本 的 标准 偏差 . 则 有 


四 | 
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如 果 样 本 数据 划分 为 m 组 ， 对 应 第 1 组 的 统计 参数 为 


e ] & 
X,2—»X, (15) 
M. i-i 


i 


2 Dy? v2 
V; => xX; -mX ’ (16) 


i-l 


S? =V? Km -1) a7 


对 应 m 组 的 总 方差 为 ，V2 Y V? ag 
i=l 
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5 Jurkevich 方法 分 析 PKS 0735+178 B 波段 光 变 周期 
Fig.5 Jurkevich method analysis the variability 


period of PKS 0735+178 in B band 
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图 6 Jurkevich 方法 分 析 PKS 0735+178 V 波段 光 变 周期 
Fig.6 Jurkevich method analysis the 


variabilityperiod of PKS 0735+178 in V band 


应 用 上 述 Jurkevich 方法 获得 PKS 0735+178 的 B 和 VV 波段 p-y, 图 ， 由 图 5、 图 6 可 


以 看 出 , 对 应 于 V 的 最 小 值 , B 波段 存在 1. 02 年 、1. 20 年 、1. 66 年 、4. 77 年 的 光 变 周期 ， 


V 波段 周期 存在 0. 94 年 、1. 71 年 、2. 82 年 、3. 97 年 的 光 变 周期 。 根据 光 变 周 期 存在 的 判 据 
m9,B 波 段 存在 1.02 年 和 VV 波段 存在 0.94 年 光 变 周期 信息 的 可 靠 性 最 大 , 而 B 波 段 存在 4. 77 


和 波段 存在 3. 97 年 次 之 . 但 考虑 到 平均 时 标 为 (1.11 土 0.01) 年 的 周期 很 可 能 是 由 于 地 球 
绕 太 阳 所 致 , 则 B 波段 存在 4. 77 年 周期 ，V 波段 存在 3. 97 年 周期 。 


2.3 离散 相关 分 析 法 分 析 


离散 相关 分 析 法 是 用 来 分 析 两 组 离散 数据 相关 性 的 方法 之 一 ”“”, 文 [24] 用 此 方法 分 析 BL 
Lac 天 体 PKS 0537-441 的 光 变 特性 ”“， 该 方法 最 大 的 特点 是 不 需要 对 数据 做 任何 处 理 ， 就 
能 判断 出 两 组 数据 的 相关 性 。 该 方法 的 具体 的 步骤 如 下 : 


首先 计算 两 组 数据 的 离散 相关 函数 值 。 如 数组 Q; 和 已 ， 则 离散 相关 函数 值 为 : 
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Ab 分 别 为 两 组 数据 的 平均 值 ，6,6, 分别 为 对 应 的 标准 偏差 。 


在 区 间 T 土 AT/2 中 有 M 个 At. 则 DCF (t ËX 


1 
DCF(t)= y LUCK, (20) 
再 次 离散 相关 函数 的 误差 为 


oDCF(t )= i [upcr, - DCF( 


1 


2d: (QD 
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其 次 计算 DCF (t ) fü. 通过 时 间 延 迟 A 方 =t,-t,; 把 两 组 数 联系 起 来 , 假如 时 间 延 迟 为 下 ， 


对 于 所 得 到 的 离散 相关 图 ， 如 果 峰 值 在 零 的 右面 ， 表 明 数 组 Q; TBD, 的 变化 。 反 之 ， 


BAA RAIL BLED, 。 


B: 1 1 1 1 1 
-1500 -1000 -500 0 500 1000 1500 


delay(day) 


delay(day) 


图 7 PKS 0735+178 在 B 波段 的 相关 性 分 析 


4] 8 PKS 0735+178 在 V 


Fig.7 The correlation analysis of PKS 


0735+178 in B band 


存在 周期 1.21 4E. 2.22 年 、3. 15 年 、4. 00 年 的 光 变 周期 信息 。 


Fig.8 The correlation 


波段 的 相关 性 分 析 


analysis of PKS 


0735+178 in V band 
通过 使 用 离散 相关 函数 对 PKS 0735+178 的 B 波段 图 像 进行 分 析 ， 
从 图 中 得 到 了 波段 周期 存在 1. 21 年 、 2.27 年 、3. 62 年 、3. 81 年 的 光 变 周期 信息 ，V 波段 


所 得 结果 如 图 7、8， 


根据 离散 相关 示 数 分 析 周 


期 的 方法 ””,B V 波段 存在 1. 21 年 光 变 周期 信息 的 可 靠 性 最 大 , M B 波段 存在 4.77 和 V 
波段 存在 3. 97 年 次 之 .但 考虑 到 平均 时 标 为 (1.11 土 0.01) 年 的 周期 很 可 能 是 由 于 地 球 绕 太 
FH Pr, JJ PKS 0735+178B 波段 的 周期 为 3. 81 FE, V 波段 的 周期 为 4. 00 年 。 


3 讨论 与 结论 


本 文 主要 收集 了 PKS 07354178 的 光 变数 据 ， 运 月 


时间 补偿 离散 傅 旦 


有 叶 法 变换 、urkevich 
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方法 和 离散 相关 分 析 法 3 种 方法 分 析 PKS 0735+178 的 光 变 周期 ， 时 间 补 偿 离 散 傅 里 叶 法 变 
换 法 分 析 B 波段 得 到 的 周期 大 致 为 5.80x10“(1/day), V 波段 周期 为 5.80x104(1/day); 


Jurkevich 方法 分 析 得 到 B 波段 的 周期 大 致 为 1740 R, V 波段 的 周期 大 致 为 1450 R; DCF 
方法 分 析 的 得 到 B 波段 的 周期 大 致 为 1390 R, V 波段 的 周期 大 致 为 1460 R. 由 此 可 得 出 ， 


BL Lac 天 体 PKS 0735+178 的 光 变 周期 为 (4.33 土 0.41)yr 。 这 与 文 [25] 和 文 [11] 用 其 他 方 


法 研究 结果 基本 一 致 。 根 据 周期 分 析 获 得 的 PKS 0735+17 的 周期 (4.33 土 0.41)yr ， 我 们 利 


用 薄 吸 积 盘 理 论 ” 来 分 析 该 赐 虎 天 体 PKS 0735+178 的 中 心 黑洞 质量 和 薄 吸 积 盘 的 热 不 稳定 
性 发 生 的 区 域 。 一 般 情况 下 ， 热 有 限 循 环 周期 性 取决 于 粘 滞 度 参数 QQ ， 中 心 黑洞 质量 为 


以 和 广义 应 力 张 量 参数 内 ， 因 此 爆发 时 间 为 
* 10°Mo 


-0. 2 M S 37 
T Durst = 4. 520, 1 YTS 


, (22) 


其 中 , OL 为 粘 滞 系 数 ; M ,以 10 ' 倍 太阳 质量 为 单位 ; M 为 PKS 0735+178 的 中 心 黑洞 质量 
le HF B Bf PARE RE DA ED RE INE, 用 磁 流体 力学 对 其 讨论 是 一 种 常用 的 
方法 , 文 [26] 提 出 如 果 磁 场 的 逃逸 率 比 较 低 ， 可 以 认为 由 = 0.5， 使 用 这 个 参数 值 ， 长 周 其 


爆发 时 间 大 约 为 2 S, B 


burst 


=> _ -0.62 4 1.37 
N bye = 9.0405 M6 YTS 023) 


对 于 PKS 0735+178， 如 果 Q = 0.1， 几 = 0.5 ,分析 获 得 的 周期 是 (4.33 土 0.41)yr， 可 


tM =0.22x10 Mo; 这 个 中 心 黑洞 质量 对 于 PKS 0735+178 来 说 


i 
= 
v 


UR HH SEU RR 


o 太 小 ， 由 于 利用 薄 吸 积 盘 理论 分 析 方 法 考虑 黑洞 的 自 旋 ， 所 以 中 心 黑洞 质量 偏 小 , 但 蝎 虎 天 
体 的 观测 特性 表明 , 这 类 天 体 在 光谱 观测 中 无 发 射线 "", 作为 深入 研究 该 类 天 体 给 出 了 黑洞 
质量 的 下 限 也 是 有 意义 的 , 在 今后 的 CCD 测 光 观测 中 寻找 短 时 标 光 变 ， 从 而 获得 中 心 黑洞 质 

量 ， 与 之 做 进一步 比较 。 

从 图 上 也 可 看 出 不 同 的 方法 周期 性 的 明显 程度 也 不 一 样 。 时 间 补 偿 离 散 传 里 叶 变 换 法 可 以 运 

到 天 体 数据 处 理 上 , 可 以 有 效 地 处 理 不 均匀 的 天 文 数据 ,克服 传统 的 傅 里 叶 变 换 带 来 的 伪 

周期 ， 时 间 补 偿 离散 傅 里 时 变换 法 通过 施 密 特 正 交 化 ， 有 效 地 解决 由 于 1、COS Qf 和 

SIN Of 3 个 向 量 不 正 交 而 带 来 伪 信 号 的 问题 ”“。 用 离散 相关 性 函数 对 PKS 0735+178 天 体 

Hj B 波段 和 7V 波段 相关 性 分 析 ， 从 图 中 可 以 看 出 它们 之 间 存 在 很 强 的 相关 性 ,表明 它们 之 间 

的 辐射 起 源 是 相似 的 ”” 。 


另外 ， 在 光 变 周 期 信息 分 析 中 ， 注 意 到 (1.11 土 0.01) 年 这 个 结果 可 能 和 一 年 的 观测 周期 有 


关 ， 这 个 周期 可 能 是 由 于 有 规律 的 观测 时 间 间 隔 导 致 的 ”， 平 均 时 标 为 (1.11 土 0.01) 年 的 


周期 很 可 能 是 由 于 地 球 绕 太阳 所 致 ”。 当然 , 也 可 能 是 天 体 PKS 0735+178 的 真实 光 变 周期 ， 
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其 他 光 变 周期 信息 是 这 个 周期 的 N 倍 (N=1，2，3…) “”。 这 个 周期 的 研究 ”"， 对 今后 的 观 
测 有 指导 意义 。 将 通过 实测 和 更 多 的 观测 证 明天 体 PKS 0735+178 的 光 变 周期 。 
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Abstract:In this paper, we used the time compensation discrete Fourier transform method and 
discrete correlation analysis (DCDFT), Jurkevich ,DCF model to analyses B band and V band 
periodicity of the PKS 0735 + 178, based on a large amount of data we have collected. We found 


that the variability period is about (4.33 EE 0.41) years and the lower limits of the central 
black hole mass is Q.22x 10* M. 


Keywords: PKS 0735 + 178; DCDFT method; Jurkevich method; DCF model method; 
Variability; Black hole mass 


